#include "mssm_xs_tools.C"

double XSCalculator(double mA, double tan, TString pwd="./")
{

//gROOT->ProcessLine(".L mssm_xs_tools.C");

mssm_xs_tools mssm_xs;
TString file = pwd + "out.mhmax_7.root";
mssm_xs.SetInput(const_cast<char*>(file.Data())); 

double result=-1;

if (mA<=120e0) result=mssm_xs.Give_Xsec_bbA4f(mA,tan)*mssm_xs.Give_BR_A_bb(mA,tan)+mssm_xs.Give_Xsec_bbh4f(mA,tan)* 
mssm_xs.Give_BR_h_bb(mA,tan);
	
else
result=mssm_xs.Give_Xsec_bbA4f(mA,tan)*mssm_xs.Give_BR_A_bb(mA,tan)+mssm_xs.Give_Xsec_bbH4f(mA,tan)*
mssm_xs.Give_BR_H_bb(mA,tan);

result/=1000;

cout<<"XS("<<mA<<","<<tan<<") ,pb= "<<result<<endl;

return result;
}
